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. ■ 1. INTRODUCTION 

> ■ 

From the observational evidence of a few strongly magnetized neutron stars, which are the sources of anomalous 
X-rays and soft gamma rays, also called magnetars UK &H, the study of crustal matter, in particular the outer 
crust of such compact stellar objects have gotten a new dimension. These exotic objects are also called anomalous 
X-ray pulsars (AXP) and soft gamma repeaters (SGR). The outer crust of a typical neutron star in general, is mainly 
p/j composed of dense crystalline metallic iron. The density of such metallic crystalline matter is ~ 10 _4 po) where p is 
' the normal nuclear density (~ 2.8 x 10 14 gm/cc). Therefore, it is absolutely impossible to investigate the properties of 
such ultra-dense matter in material science laboratories, even at zero magnetic field. The observed surface magnetic 
field of the magnetars is ~ 10 15 G, which is again too high to achieve in the terrestrial laboratories. Also, it is quite 
I , possible that the interior field of such exotic objects can go up to ~ 10 18 G (which can be shown theoretically by Virial 
theorem). If the magnetic field at the interior is really so high, then most of the physical and chemical properties 
of the dense neutron matter should change significantly from the conventional neutron star (radio pulsars) scenario 
(see the recent article by one of the co-authors of this article [5[ for necessary references). Even though the magnetic 
i field strength at the crustal region is slightly higher than 10 15 G, it must change significantly most of the properties 
of dense matter, both in the outer crust and the inner crust regions of the magnetars 0, Q. It is believed that strong 
magnetic field can cause a structural deformation of the metallic atoms present in the outer crust of a neutron star. 
The spherical symmetry of the atoms will be destroyed and becomes cigar shape with the elongated axis along the 
direction of strong magnetic field. The atoms may even become almost an one dimensional string like object, i.e., 
needle shape, if the magnetic field strength is extremely strong. In a future article we shall present the problem related 
y—{ • to structural deformation of atoms in a strong quantizing magnetic field using a completely different approach 0]. 

In this article we shall study the properties of outer crust matter composed of magnetically deformed metallic (iron) 
atoms. In section 2 we have developed the basic formalism and discuss the numerical results, whereas in the last 
' section we have given the conclusion and the future prospect of this work. In this article, for the sake of simplicity, we 
shall assume a cylindrical deformation of the atoms in the outer crust region and use cylindrical coordinate system 
with azimuthal symmetry. In reality, to investigate the cigar like deformed atoms in presence of strong magnetic field 
one has to use prolate spheroidal coordinate system. 
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BASIC FORMALISM 



The width of the outer crust of a typical neutron star is ~ 0.3 km, the density of matter, which is assumed to be a 
dense crystal of metallic iron is ~ 4.3 x 10 1:L gm/cc @. To investigate the properties of such dense exotic crystalline 
matter of metallic iron, we have replaced the outer crust matter by a regular array of cylindrically deformed (in our 
future work we shall consider more realistic geometrical structure of the deformed atoms, which is prolate spheroidal 
in nature) Wigner-Seitz (WS) cells, with the positively charged nucleus at the centre surrounded by a non-uniform 
electron gas. The axis of each cylinder is along the direction of magnetic field and further assumed azimuthal symmetry 
for all the cylindrically deformed WS cells. 

We start with the Poisson's equation, given by 

V 2 = 47ren e (1) 

where 4> is the electrostatic field, e is the electron charge and n e is the electron density, which because of assumed 
non-uniformity within the WS cell, is a function of positional coordinates (r, z) . Now in the cylindrical coordinate 
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with circular symmetry, the above equation reduces to 

d 2 (t> ld<j> d 2 4> 

ttt + - tt + ttt = 47ren e (2) 
or z r or az l 

It is well known that in the presence of strong quantizing magnetic field, the number density of degenerate electron 
gas is given by 

eB Vmax 

i/=0 

where B is the constant external magnetic field, assumed to be acting along Z-direction and is > B^\ where is 
the typical strength of magnetic field beyond which, in the relativistic region the Landau levels for the electrons are 
populated. For the sake of convenience, throughout this article we shall use % = c = 1. The critical strength is given 
by B c = rr? e j | e | Q and assume that the matter is at zero temperature, where m e is the electron rest mass and 
| e | is the magnitude of electron charge. In eqn.(3), pp is the electron Fermi momentum, v is the Landau quantum 
number, with v m ax, the upper limit of v. The factor (2 — S v o) takes care of singly degenerate v = state and doubly 
degenerate all other states with v ^ 0. To study the properties of outer crust matter with deformed WS cells, we make 
Thomas-Fermi approximation, which is a semi-classical approach and having a lot of shortcomings and limitations 
@ . We further use the Thomas- Fermi condition, given by 

Me = (p F + m l + ZveB) 1 / 2 — ecf> = constant (4) 

where fj, e is the electron chemical potential and in this model, it is assumed to be constant throughout the WS cell. 
From this equation we can express the Fermi momentum of electrons in the following form: 

PF = [([i e + ec/)) 2 - m 2 e - 2veB] 1/2 (5) 

Since the electrostatic potential <j> = tp(r,z), the Fermi momentum pp for electron is also a function of positional 
coordinates (r, z) within the cell. In principle one should use the exact expression for electron Fermi momentum as 
given above, in the equation for electron density (eqn.(3)) which is in turn appearing on the right hand side of the 
cylindrical form of Poisson's equation (eqn.(2)). However, with this exact expression for pp, it is absolutely impossible 
to proceed further analytically, even a single effective step. From the very beginning, therefore, one has to use some 
numerical technique to solve the Poisson's equation, which is an extremely complicated nonlinear partial differential 
equation. Of course, with the numerical method within the limitation of the algorithm followed, we will get more 
exact results. However, in numerical computation of cf>(r,z), we get a set of numbers, but the beauty of this model 
will bc completely destroyed and a lot of interesting physics associated with intermediate results of this problem will 
be totally lost. In this context, let us quote from a quite old but extremely interesting paper by Wigner and Bardeen 
[Tfj| - "It is perhaps not quite superfluous to have, in addition to a more exact calculations of a physical quantity, an 
approximate treatment which merely shows how the quantity in question is determined, and the lines along which 
a more exact calculation could be carried out. Such a treatment often leads to a simple formula by means of which 
the magnitude of the quantity may be readily determined" . Therefore, to get an approximate analytical solution for 
ip(r, z), we assume that v max = and neglect the rest mass of electron, i.e., we put m e — in the expression for 
Fermi momentum pp (eqn.(5)). The approximation v max — is actually valid if the magnetic field is extremely high. 
However, to investigate the properties of dense electron gas within the cylindrically deformed cells in a little bit of 
exact manner later in this article, we shall use the approximate solution for cf>(r, z), but do not restrict ourselves to 
Vmax = 0. In fact we shall show later that the upper limit v max for the electron Landau quantum number is also a 
function of (r, z) . Therefore to evaluate various physical quantities in the outer crust region, we have to first obtain 
an approximate solution for the Poisson's equation. To achieve our goal, we use the approximate form of electron 
Fermi momentum obtained from the assumption as mentioned above and its mathematical form is given by 

p F ~ n e + ecj> (6) 

Next on substituting 

H e + e<t>{r,z)=ij){r,z), (7) 
the cylindrical form of Poisson's equation reduces to 

d 2 ip 1 dtp d 2 i/j 



dr 2 r dr dz 



, - X ip (8) 
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where A 2 = 2e 3 B/ir. Under this approximation, the Poisson's equation, as shown above, reduces to a linear partial 
differential equation. To solve this partial differential equation analytically we use the method of separation of 
variables, given by 



ip{r,z) = R(r)Z(z) 

Substituting z) from eqn.(9) in eqn.(8) and introducing a constant £, we get 



CR = 



d 2 R 
dr 2 



1 dR 

r dr 

f 2 , >2 



(9) 

(10) 
(11) 



where £ is some real constant, independent of r and z but may change with the magnetic field strength and with 
the mass number and the atomic number of the type of elements present in the outer crust region. The solutions of 
eqns.(lO) and (11) are well known. For eqn.(10), the solution is a Bessel function of order zero with the argument £r, 
whereas for eqn.(ll), it is an exponentially decaying function of z. In the language of mathematics, the solution for 
ip(r, z) is then given by 



V>(r, z) = CJ (^r)exp [±(£ 2 + \ 2 fl 2 z 



(12) 



where + and — signs are for z < or > respectively. Here C is a constant (again may change with the magnetic 
field strength and with the atomic properties of the elements present in the outer crust) and Jnjjr) is the Bessel 
function of order zero. Now on the nuclear surface, at the centre of the WS cells, <f> = Zoe/r n [11(, where Zq is the 
atomic number and r n = tqA 1 ^ is the nuclear radius, tq — 1.12fm and A is the mass number (here we have assumed 
that the nuclei remain spherical, even if the magnetic field is too high). For the sake of simplicity we put z = and 
r = r n on the nuclear surface, i.e., we have chosen points on nuclear surface along r-axis, which gives 



ip(r n ,0) = CJ Q (£r n ) = fi e + 



Z e 2 



(13) 



Hence we get 



C = 



Jo(£r n ) 



fj, e 



Z e 2 



(14) 



In the expression, of course, £ is still an unknown quantity. To determine £, we now put z = +r n and r — in the 
expression for tf>(r,z) on the nuclear surface, i.e., we have chosen a point on the positive side of Z-axis. Then using 
Jo(0) = 1, we have (because of the symmetry about z — plane, z — ~r n and r — will also give the same result) 



^(0,r„) = Cexp[-(£ 2 + A 2 ) 1/2 r, 



= He 



Z e 2 



(15) 



Hence we get 



exp 



. ( £2 + a2) 1/2 1 =Jo ^ rn) 



(16) 



This is a highly transcendental equation for £. However, it is possible to evaluate £ numerically from this equation 
for a given magnetic field strength and for a given type of element, say metallic iron in the crystalline form present 
in the outer crust region. The above expression puts some constraint on the range of £r„. To make this point more 
transparent let us express eqn.(16) in the following form 



£ + A 



1 



[In J (£r„ 



(17) 



Since Jo(£r n ) is an oscillatory function and becomes negative as £r n crosses the first zero, in the numerical computation 

we therefore put a check in the computer code to restrict £r„ from 0+ to the first zero—, which is 2.40482 In fig.(l) 

we have plotted £ (in MeV) as a function of magnetic field strength B, expressed in terms of critical magnetic field 

strength . The variation is insensitive for the low and moderate values of magnetic field strengths. However, the 
overall magnitude of £ is very close to lOOMeV. This figure shows that beyond field strength 10 17 G, £ increases sharply. 
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In fig. (2) we have plotted the variation of normalization constant C with the magnetic field strength, expressed in 
the same unit as in fig.(l). For low and moderate magnetic field values it is almost constant, then it falls abruptly 
beyond 10 17 G and finally saturates to a constant value equal to ~ lOOMeV. Since beyond 10 17 G, electrons within the 
cells occupy their zeroth or very low lying Landau levels, the quantum mechanical effect of magnetic field becomes 
extremely important and as a consequence both £ and C change significantly beyond this magnetic field value. 

Since the WS cells are overall charge neutral, then at any point (r, z) on the surface, even on the plane faces or on 
the boundary points of z = plane, we have in cylindrical coordinate 

= (18) 

This condition gives 

e;?+4?=0 (19) 



dr d 



z 

which must be satisfied at all the surface points on the WS cylinder. 
Now from the Thomas-Fermi condition, we have 

p F = ^ 2 {r,z)-m 2 ] 1 / 2 (20) 

where m„ = (to 2 + 2veB) 1 / 2 . With this more exact expression for pp, the number density for electron gas can be 
expressed as 

n e (r, z) = g £ (2 - M[V>V, z) - ml] 1 ' 2 (21) 

This is obviously more exact than eqn.(3). Further, this expression shows that the electron density is a function of 
both r and z within the WS cell. Which justifies the assumption that the electron distribution inside each WS cell 
around the fixed nucleus is non-uniform. Now from the non-negative nature of we have 

_ i> 2 (r,z)-m 2 e 

^max — 2eB ~ max ^ ' \ I 

The upper limit of Landau quantum number v max will therefore also depend on the positional coordinates of the 
associated electron within the cell. In fig. (3) we have plotted electron number density in terms of normal nuclear 
density multiplied by 10 4 , as a function of radial distance from nuclear surface to the WS cell boundary for the 
magnetic field strengths 10 1 , 5 x 10 2 , 5 x 10 3 , 10 4 and 5 x 10 4 times B c e \ indicated by the curves a, 6, c, d and e 
respectively. These curves show that the electron number density is maximum near the nuclear surface and minimum 
near the WS cell boundary r max . This figure also shows that the value of electron number density increases with the 
strength of magnetic field. In fig. (4) we have plotted the same quantity as in fig. (3) but against the axial distance 
z from the nuclear surface to the WS boundary indicated by z max . In this case also the variations are exactly same 
as in fig. (3). The qualitative difference is because of different types of functional dependence. In fig. (5) we have 
plotted the upper limit of Landau quantum number v max as a function of radial coordinate from the nuclear surface 
to the WS cell boundary. In this figure the upper curve is for 10 x Bc C \ the middle one is for 10 2 x B^ and the 
lower one is for 500 x Bi e \ The value of v max decreases with the increase in magnetic field strengths. It has been 
observed that beyond 500 x B^ , v. max becomes identically zero throughout the WS cell. Further, the value of v. max 
for B < 500 x B { c e) is largest near the nuclear surface and exactly zero at the cell boundary. In fig. (6) we have plotted 
the same kind of variation of i^rnax-, 

but against the axial distance from nuclear surface to the cell boundary. The 
qualitative and the quantitative variations are exactly identical with fig. (5). These two figures show that the electrons 
are completely spin polarized in the direction opposite to the magnetic field B at the cell boundary for B < 500 x Be , 
but beyond this value they are polarized at all the points within the cell. Although the variations along radial and 
axial directions are shown in these two figures, we expect that such polarized picture of electron gas will be there 
throughout the cylindrically deformed WS cell surface, including the two plane faces. 
Now the minimum possible value of v max is zero, hence we have 

ip 2 (r, z) — m 2 > 

To get the length and radius of the cylindrically deformed WS cell, we first put r = in the above expression. Then 
we have 



C 2 exp[±2z(A 2 +£ 2 )] > m 2 



(23) 
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Therefore beyond some value of z, say z max , this inequality will break, i.e., just at value z = z max , we have 

C 2 cM±^ma X ^ 2 +e) 1/2 }=m 2 e (24) 
Since z-axis is symmetric about z = plane, we have from this relation 

' -Inf^l (25) 



max (A 2 + C 2 ) 1/2 VC 

This is the magnitude of z on the plane faces of a WS cell, i.e., 2z max is the length of the cylindrical^/ deformed WS 
cell for a given B and a set of (Z ,A). In fig. (7) we have shown the variation of z max with the strength of magnetic 
field B expressed as before in terms of critical field strength This figure shows that the variation is almost 

insensitive for low and moderate magnetic field strengths but decreases almost abruptly beyond 10 16 G, when most of 
the electrons occupy their zeroth Landau level, and finally tends to saturate to a constant value ~ lOfm. 
To get the radius of such cells, we next put z = 0, then we have 

V> 2 (r,0) - ml > 

which further gives for r = r max , the radius of the cylindrical WS cells, 

77% 

Jo(£r max ) = (26) 

Obviously an analytical solution for r max is absolutely impossible, we obtain numerically the value of r max for a given 
B and (A,Z ). In this case also, as before, we restrict the value of the product £r max < the first zero of Jo(£,r max ). 
The variation of r max with the strength of magnetic field is shown in fig. (8). The nature of variation is more or 
less same as that of z max - However, we have noticed that for extremely large field strength, r max — ► 0. This is 
a remarkable difference from its longitudinal counter part. It actually shows that in presence of extremely strong 
magnetic field the cylinders become more and more thin in the transverse direction. We therefore conclude that with 
the increase in magnetic field strength the radial contraction will be enormous compared to the axial one. From 
figs.(7)-(8) we have noticed that the variations are most significant beyond B = 10 16 G. The reason is again because of 
the fact that the electrons occupy only their zeroth Landau level in presence of such strong magnetic field, at which 
the quantum mechanical effect of the magnetic field dominates. 

Next we calculate the different kinds of energies of electron gas within the WS cells. The cell averaged kinetic 
energy density of an electron is given by 

/ dpM+™>l) 1/2 -me] (27) 

Jo 



where in the cylindrical coordinate system with azimuthal symmetry d 3 r = 2irrdrdz, with the limits r n < r < r max 
and r n < z < z max and V = Trrf nax 2 x z max , the volume of each cell. The p z integral is trivial and can be obtained 
an analytic expression for the average kinetic energy density. However, if we do not integrate over r and z, then from 
this equation we get the energy density £fe(r, z) at a particular point inside the cell. Of course, in that case we need 
not have to divide by V. The factor 2 is for z-symmetry about z — plane. In fig. (9) we have plotted the variation of 
kinetic energy density as a function of radial distance r from the nuclear surface to the cell boundary, keeping z = 0. 
This figure shows that the kinetic energy density increases with the increase in magnetic field strength. We have 

indicated the curves by a, b, c, d and e for the magnetic field strengths \QBc \ 5 x 10 2 Bc £ \ 5 x 10 3 £>c 6 \ lO^B^ and 
5 x 10 4 i?c respectively. In fig. (10) we have shown the same kind of variation along axial distance from the nuclear 
surface to one of the plane faces of the WS cell. The variation with magnetic field strength is again almost identical 
with fig. (9). Similar to the variation of electron number density within the cell (fig. (3) and fig. (4)), both fig. (9) and 
fig. (10) show that the kinetic energy density for electrons is maximum near the nuclear surface and minimum at the 
cell boundary. 

Similarly, the cell averaged electron-nucleus interaction energy per unit volume is given by 

E ---V Zoe2 J d3r VT^ (28) 

Analogous to the kinetic energy, here also one can obtain the interaction energy per unit volume, E en (r,z) at a 
particular point within the cell. In fig. (11), the variation of the magnitude of electron- nucleus interaction energy per 
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unit volume with the radial distance is plotted. In this figure we have indicated the curves by a, 6, c, d and e for the 
magnetic field strengths lOB { c e) , 5 x 10 2 ^ e) , 5 x 10 3 B^ e) , 10 4 B^ e) and 5 x 10 4 M e) respectively. This fi gure shows that 
the magnitude of electron-nucleus interaction energy increases with the increase in magnetic field strength. Which 
means that with the increase in magnetic field strength the electrons become more strongly bound by the nuclear 
Coulomb attractive potential. It is also obvious from the figs.(9)-(12) that for a given magnetic field strength the 
kinetic energy density and the magnitude of electron- nucleus interaction energy per unit volume at a particular point, 
either along axial direction or in the radial direction, inside the WS cell are of the same order of magnitude. We do 
believe that this is true at all the points inside the WS cell. In fig. (12) we have shown the same kind of variation 
along z-axis. Both the qualitative and the quantitative nature of variations with the strength of magnetic field are 
same as that of fig. (11). 

Next we consider the cell averaged electron-electron direct interaction energy density, given by 

Efl = ye 2 J d*rn e (r,z) J rfVn e (r', z') [(f _ p)2 + ^ _ ^ )2]1/2 (29) 

Now assuming r as the reference axis, we have (r — r 1 ) 2 = r 2 + r' 2 — Irr' cos#, where 6 is the angle between r and 
r*. Then d 3 r' = r'dr'dOdz' , with < 8 < 2ir. In this case the z integral has to be broken into two parts, one with 
limit — Zmax < z < z and the other one with the limit z < z < +z max . The value of z is not easy to evaluate 
in the region between Z-axis and r-axis. With 9 symmetry, for the sake of simplicity we put z — r n and expect 
that the error will be nominal. Now to obtain electron-electron direct interaction energy one has to evaluate the five 
dimensional integral as shown in eqn.(29). None of them can be obtained analytically, hence it is necessary to follow 
some numerical method. Even the 9 integral can not be obtained analytically. One can express the 9 integral in the 
form of an elliptical integral of first kind, given by 

hc{ry,^) = ['\e {i _ K l os2e)l/2 (30) 

where K = Arr' /[(r + r') 2 + (z — z') 2 ]. The direct part is then given by 

Efl = W*f d " rn ^ J r'dr'dzn e (r>, z>) [(f + ^ + ^ _ ^ )a]1/a he (r, r', z, z') (31) 

where the factor 4 is coming from the angular integral over 6 from to 27r. Let us now consider the elliptic integral 
(eqn.(30)) onz = 0-plane. In this case both z and z' are zero and the factor K = 1. Then it can be shown very 
easily that the integral given by eqn.(30) will diverge at the lower limit. To avoid this unphysical infinity we put a 
lower cut off (infrared cut off) 8, which will now the lower limit for 9. The physical meaning of non-zero lower limit 
for the ^-integral is that the two electrons under consideration can not be at zero distance from each other on the arc 
of a circle whose centre is same as that of the nucleus. The infrared cut off 6 which is a measure of angular distance 
between two neighboring electrons must necessarily depends on the minimum possible linear distance between them 
and also on the radial distance from the centre. From a very elementary geometrical construction it can be shown 
that 

S = - (32) 
r 

where s is the arc length, or the distance between two neighboring electrons on the circular arc. Since s is infinitesimal 
in nature, we can approximate it by a straight line of length s, which is the length of the cord connecting two points 
occupied by two neighboring electrons. Since s is the minimum possible distance between two electrons, we can 
express it in terms of electron density near those points, given by 



7? 



1/3 



(33) 



In a future publication we shall evaluate this multidimensional integrals along with the electron-electron exchange 
energy (as given below) using the spinor solutions of Dirac equation in cylindrical coordinate in presence of strong 
quantizing magnetic field using a sophisticated Monte-Carlo numerical code for multidimensional integrals [1 21 ] - 
The electron-electron exchange energy corresponding to the ith electron in the cell is given by 

E it x) = -y E / M \ { ,_ ??+ \ Z _ Z , YV/ J ^ *0fc(* -0 (34) 
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where ipi{r, z) is the spinor wave function in cylindrical coordinate in presence of strong quantizing magnetic field, and 
z) — rji iT(r, z)j , the adjoint of the spinor and 70 is the zeroth part of the Dirac gamma matrices 7^ in cylindrical 
coordinate system. 

The kinetic pressure of non-uniform electron gas within the WS cell is given by 

P(r,*) = |4 £ (2-M/ , 2 P r% /2 (35) 

The p z integral is very easy to evaluate analytically and is given by 

u=0 



- m llJ PF + ^ + mv2)1/2 \] (36) 

This equation shows that the electron kinetic pressure also changes from point to point within the WS cells. In fig. (13) 
we have shown the variation of kinetic pressure for the non-uniform electron gas with the radial distance within the 
cell. Curves a and b are for B = 10 x and B — 500 x B^\ whereas upper curve and lower curve as indicated by c 
(almost identical) are for B = 5000 x B^ and B = 10000 x B^ respectively. In fig. (14) the same kind of variations 
are shown along z-axis. In this figure the curves for B = 5000 x B [ c ] and B = 10000 x B ( c e) are almost identical 
and indicated by single thick curve c. For B = 10 x B [ c ] and B = 500 x Bi e) the curves are indicated by a and b 
respectively. These two figures show that the kinetic pressure is maximum near the nuclear surface and zero at the 
cell boundary. The variation with magnetic field strength shows that the non-uniform electron gas becomes softer for 
high magnetic field. 



3. CONCLUSIONS 



In this article we have investigated various physical properties of non-uniform electron gas within the cylindrically 
deformed atoms of metallic iron at the outer crust of a strongly magnetized neutron star (magnetar). Because of 
extremely strong magnetic field, we have assumed a cylindrical type deformation of the atoms, replaced by WS cells 
with the same kind of deformation. The axes of the cylinders are along the direction of magnetic lines of forces. The 
curved surfaces of these cylinders are therefore parallel to the boundary surface of the neutron stars in the region 
far away from the magnetic poles. We have studied the variation of electron number density, kinetic energy density, 
electron-nucleus interaction energy per unit volume of electron gas and also kinetic pressure of non-uniform electron 
gas within the cylindrically deformed WS cells with the strength of magnetic field for a given spatial coordinate (r, z) 
and also with the radial and axial distances within the cell for a given magnetic field value. We have also investigated 
the variations of longitudinal and transverse dimensions of deformed WS cells with the strength of magnetic field. We 
have noticed that the transverse dimension of a cylindrically deformed WS cell becomes extremely thin in presence 
of ultra-strong magnetic field. It is believed that the atoms become cigar shape in presence of strong magnetic field. 
In our future study, the properties of neutron star crustal matter with cigar shape atoms in the metallic crystal in 
presence of strong magnetic field. In future, we shall also evaluate the electron-electron direct interaction energy and 
the exchange part of electron-electron interaction. 
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FIG. 1: The variation of separation variable £ with the magnetic field strength 




FIG. 3: The variation of electron number density expressed in terms of normal nuclear density with radial distance r in Fermi. 
The curves indicated by a, b, c, d and e are for B = 10, 5 x 10 2 , 5 x 10 3 , 1 x 10 4 and 5 x 10 4 times B< e) respectively. 
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FIG. 4: The variation of electron number density expressed in terms of normal nuclear density with radial distance z in Fermi. 
The curves indicated by a, b, c, d and e are for B = 10, 5 x 10 2 , 5 x 10 3 , 1 x 10 4 and 5 x 10 4 times B ( c e) respectively. 
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FIG. 5: The variation of upper limit of Landau quantum number with the radial distance in Fermi. Upper curve is for 
B = 10 x B ( c e) , middle one is for 10 2 x B ( c e) and the lower one is for 5 x 10 2 B^ e) . 




FIG. 7: The variation of z max with the magnetic field strength. 
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FIG. 10: The variation of kinetic energy density with the axial coordinate in Fermi. The curves indicated by a, b, c, d and e 
are for B = 10, 5 x 10 2 , 5 x 10 3 , 1 x 10 4 and 5 x 10 4 times B ( c e) respectively. 
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FIG. 11: The variation of electron nucleus interaction energy density with the radial distance in Fermi. The curves indicated 
by a, b, c, d and e are for B = 10, 5 x 10 2 , 5 x 10 a , 1 x 10 4 and 5 x 10 4 times B ( c e> respectively. 
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FIG. 12: The variation of electron nucleus interaction energy density with the axial coordinate in Fermi. The curves indicated 
by a, b, c, d and e are for B — 10, 5 x 10 2 , 5 x 10 3 , 1 x 10 4 and 5 x 10 4 times Bi e) respectively. 
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FIG. 13: The variation of kinetic pressure of non-uniform electron gas within WS cell with the radial distance in Fermi. Curves 
a and b are for B = 10 x Bi^ and 5 x 10 2 Bc e ' respectively, whereas upper and the lower curves indicated by c are for 5 x 10 3 B^ 
and 10 4 x B^ respectively. 
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FIG. 14: The variation of kinetic pressure of non-uniform electron gas within WS cell with the axial coordinate in Fermi. 
Curves a and b are for B = 10 x and 5 x 10 2 Bc^ respectively, whereas upper and the lower curves indicated by c are for 
5 x 10 3 B^ e) and 10 4 x B ( c e) respectively. 



